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Abstract 

In the present paper, we study the non-Gaussianity of multiple-field infla- 
tion model using the method of the field equation. We start from reviewing the 
background and the perturbation theory of multiple-field inflation, and then 
derive the Klein-Gorden equation for the perturbations at second order. Af- 
terward, we calculate the tree-level bispectrum of the fields' perturbations and 
finally give the corresponding parameter f^L for the curvature perturbation 
in virtue of the 5N formalism. We also compare our result with the one already 
obtained from the Lagrangian formalism, and find they are consistent. This 
work may help us understand perturbation theory of inflation more deeply. 
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1 Introduction 



It is suggested that our universe has undergone an inflationary stage in the early 
time. This scenario helps us understand why our universe are so flat and isotropic, 
and also provides a possible solution to the monopole problem in the hot Big-Bang 
cosmology [U [2[ |3]. The most efficient model of inflation is driven by a single scalar 
field which rolls down along its potential very slowly. This model generically predicts a 
scale-invariant powr spectrum and so is able to explain the formation of the large scale 
structure. This expectation has already been confirmed by the 5-year WMAP data ^ 
which is the latest observation of cosmic microwave background radiation (CMBR). 
Although the single field inflation model has obtained fruitful achievements, we still 
need to explore more on this theory. 

A significant lesson is to investigate its higher order perturbations. There has been 
a number of literature studying the behavior of higher order perturbations in inflation 
models (see [5] for an excellent pioneer work, and see [5] for a good review on this 
issue). From the viewpoint of statistic dynamics, these higher order perturbations 
are usually related to n (n > 2) point correlators. So if these correlators indeed 
exist, there must be non-Gaussianity in the early universe. The non-Gaussianity is 
a very important issue worth studing. Since the non-Gaussianity has many features 
which can be observed by experiments, such as its magnitude, shape, running and so 
on, it encodes plentiful information about the early universe. We are able to learn 
what has happened since that time if we detect it. For example, we have already 
known that the primordial bispectrum of single scalar field inflation model is too 
small to be observed [5j. However, there are implications of non-Gaussianity which 
value may be large from astronomical data [H [7] recently. If this is confirmed, the 
usual inflation models, especially chaotic inflation, will suffer a great challenge from 
the experiments. Another example is, that different models usually have different 
predictions about non-Gaussianity [6l[Hl[9l[Tni[ni[l2l[l3l[Tll[T5l[l6j, and thus the 
non-Gaussianity may help us to discriminate these models. 

One may be interested in how to produce large non-Gaussianity in inflation. 
Firstly, let us think about why the non-Gaussianity of single scalar field inflation 
model is so small. Ordinarily, the non-Gaussian effect comes from the interactions 
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of the perturbation variables in the canonical case that we choose the Bunch-Davis 
vacuum. In the single scalar field infiation model, the interaction terms of the per- 
turbation variables are strongly suppressed by slow roll parameters. However, if we 
modify the lagrangian of the infiaton to be non-canonical, such as DBI infiation [9l |14] 
and K-infiation [10], it is possible to obtain a large value of non-Gaussianity. Another 
way to obtain large bispectrum is to change the initial condition of Gaussian statistic, 
e.g. a thermal initial condition [15] . 

The above arguments are valid when we only consider the adiabatic perturbations. 
It is feasible since the perturbations generated in single scalar field infiation model 
are always highly adiabatic and the curvature perturbation ( is conserved on large 
scale. However, this picture is changed when the infiation is driven by multiple 
fields [171 [18] . When we introduce multiple fields, they will generate a large amount 
of entropy fiuctuations which are converted into the curvature perturbations at later 
time. This scenario can result in large non-Gaussianity of local form. For example, 
the curvaton mechanism [191 120] and the in-homogenous reheating scenario [211 1221 [23] 
are able to produce large non-Gaussianity of local form. Therefore, it is meaningful 
to study the bispectrum of multiple-field infiation with both the magnitude and the 
shape in detail. 

Interestingly, the method of calculating the primordial non-Gaussianity is not 
unique, and a number of methods have been proposed in literature. These methods 
possess different advantages in different occasions. The most direct formalism was 
developed by [23] in which the authors calculated the second order perturbations from 
the Einstein equations; another useful formalism was called Lagrangian formalism [5] 
which derived the interaction terms of curvature perturbations in the Lagrangian. 
Both the two approaches are able to calculate the magnitude and the shape of the non- 
Gaussianity of which the local one is the most interested in observations. Moreover, a 
so-called 6N formalism [25] has been proposed to calculate the local non-Gaussianity 
specifically. This method greatly simplified the calculation of non-Gaussianity. Some 
pioneer works on the non-Gaussianity of multiple-field infiation have been done by 
using different methods. For example, the local form of non-Gaussianity in two-field 
infiation is shown by [26] based on 6N formalism; the shape of non-Gaussianity in 
multiple-field infiation is given by the Lagrangian formalism in [27]. Recently, a 
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remarkable work has been done by [28] in which the authors have used second-order 
Klein-Gordon equation [221 ED] to calculate the non-Gaussianity, which is consistent 
with the Lagrangian formalism. The method of field equation can directly derive 
the non-Gaussianity from equation of motion, without assuming an effective action 
principle. We in this paper extend the field equation formalism to the multiple-field 
inflation and calculate the non-Gaussianity. In the derivation, we assume that there 
are M scalar fields ip^ , (p'^ , ■ ■ ■ in the period of inflation, and the potential V of the 
scalar fields depends on them. We take the natural unit Mp = (SttG)"^/^ = 1 in this 
paper. 

Our paper is organized as follows. In Section ^ we review the background evo- 
lution of the multiple-field inflation, and define the slow-roll parameters. In section 
^ the quantum theory of the first order perturbations is discussed. By means of 
the canonical method, we quantize the perturbations of scalar fields, and present the 
Green's functions. In Section ^ we derive the second order Klein-Gordon equation 
directly from the action, and so the second order perturbations of scalar fields are 
obtained by the Green's function. Section ^ presents the main result of our paper 
which shows that there are different source terms contributing to the three-point 
correlator of scalar fluctuations. In Section ^ we review the 6N formalism, and 
calculate the nonhnear parameter /nl- Conclusions and discussions are summarized 
in the last section. 



2 The background in multiple-field infiation 

In this section, we show the field equations in the background, and define some slow 
roll parameters in multiple-field inflation. The background is assumed to be the 
Friedmann-Robertson- Walker (FRW) spacetime, and the action takes the form 

ds^ = -dt^ + a{tf5ij dx' dx^ , (1) 

where a{t) is the scale factor. In some cases, it is convenient to use conformal time 
?7, which is defined as r/ = J^°° dt' /a{t'). And to the leading order of slow roll approx- 
imation, ~ — ^ in the period of inflation. 
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The equation of scalar field takes the form 

cj>i" + 2ncl>'o' + Vi = Q, (2) 

where / denotes different scalar fields, prime denotes Vj is the shorthand for 

H = a'/a is the conformal Hubble scale, and the metric of the field space is assumed 

to be 6ij. 

The 0-0 component of Einstein equations gives the so-called Friedmann equation, 

Sn^^^Sij<Pi'<p'o' + a'V{<Po) . (3) 
And from the i-j component of Einstein equations, we have 

+ 2n' = -UiAl'(t>i' + aV(0o) , (4) 

where the repeated up index and down index represent summation. 

As in single field inflation, the potential should satisfy the slow roll condition due 
to the constraint from the observation. It requires that the velocity and acceleration 
of inflaton rolling down the potential are very small. In the multiple-fleld inflation, 
we use the slow roll matrix 

TJ _ 0000 _ 4>oK _ I J /.N 

where 

e'-^. (6) 

and the trace of the slow roll matrix tre^"' is the standard slow roll parameter e = 
—H/H^. In general situation, the order of these slow roll parameters are estimated 
as 

e"~0(i^),e'~0(^). (7) 

To generalize the single field infiation, we introduce the second slow roll matrix, 

IJ ^ 0V^ + 0^0-^ .g. 

AHH 

The diagonal element of this matrix is the slow parameter in the single field inflation 

■n'P't' — ^ = ri 

I H<j> '■ 
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3 The First-order perturbation in the uniform cur- 
vature gauge 



When we compute the perturbation of inflation, the quantity is usually changed with 
the coordinate transformation. In order to discuss the real physical freedoms in the 
inflationary perturbation theory, we should select a gauge jSTj. Fixing a gauge means 
choosing a coordinate system. Different gauges are equivalent in physics. In this 
section we select the uniform curvature gauge and discuss the first order perturbation 
of real physical freedoms. It is convenient to study in ADM formalism and the metric 
can be expressed as 

d = -NMt^ + hij (dx* + Kdt) {dx^ + Wdt) , (9) 

so the action is 

S = -Ij Nv^ {5ijh'^d4'd,<p' - 2V{<P)) 

+i j N-'Vh I^EijE'^ -E^ + - Wd,^')i^' - Wd,^')) , (10) 

where N~^Eij is the extrinsic curvature, E = E^. We select the uniform curvature 
gauge, in which the Ricci curvature is zero at the same coordinate t and hij = a?{t)5ij. 
The two scalar perturbations from the metric perturbation can be expressed by the 
lapse N, and shift N*. The lapse N, and shift N* are Lagrangian multipliers. Thus the 
physical freedoms can be expressed by the M scalar perturbations 50^ in the uniform 
curvature gauge. 

As in the single field, the scalar perturbation can be expanded in powers of the 
gaussian perturbation 

5ct>' = 5ct>{ + U<Pi + --- + ^^5ct>l + --- . (11) 

The closer the primordial scalar perturbation is to gaussian statistics, the better the 
expansion is. 

Since obeys the gaussian statistics, the equation of motion of is linear. 
After some simplification of ffTOl) . the second order action takes the form 

S2 = \j dT^d^x {5iM{'5<Pi' - SjjdScp[d5<p() , (12) 
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where dS(f){dS(l)'( is the shorthand for the scalar product 6^^di6(f){dj6(f)(. Then the field 
equation of scalar field 6(f)^ for the Fourier mode is 



i" + 2m<p{' + k'^5<i){ = 0, (13) 
The classical field is quantized by the canonical method, 

S4>i{^,v) = J ^^s^'-'-WMiv) + aLAm , (14) 

where dl, 6l are massless scalar fields in momentum space. The normalization of the 
terms is determined by the commutative relation between scalar field and its canonical 
momentum, and the commutative relation between the creation and annihilation 
operator 

[ai,ai}] = {2n)'6'-'6{k-k') . (15) 
In the Bunch-Davies vacuum, the normalized scalar field is [32], 

ei = -^{l-tkr^)e^^^ . (16) 

Since the value of is independent of /, we omit the index / in afterwards. The 
two-point correlator of scalar fields is 

(50((k,r7)505f(k',V)) = {27^f5''5{k + \^%{v)0M) 

~ {27Tf6'^6{k + k')— for kr]<^l 

= (27r)35"5(k + k')^P(fc) , (17) 

where P(/c) = ^ is the so-called power spectrum of scalar field. And the retarded 
Green's function in momentum space takes the form. 



Grk{v,r) = iairy x 



1] < T 

ek{r)h{v)-Ok{r)ek{v) V>r 



Using the Green's function, the second order field equation can be solved as a linear 
function with the source term. 
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4 The second-order Klein-Gordon equation 



In this section, we derive the second-order Klein-Gordon equation from the multiple- 
field action (fTOl) . The situation of single field is given by |30J. Expanding the action 
(ITU]) , it includes the terms of all the scalar fields and scalar perturbations from the 
metric. The lapse and the shift in the action are determined since they are Lagrangian 
multipliers without dynamics effect. Finally they are eliminated from the action which 
only contains the second-order perturbation of scalar fields 502- The part of the action 
quadratic in can be expressed in conformal time, 

S2 = ^j dr]d^x {5ij5<p{' 5<pi' - 5ijd5(t){d5<pi) . (19) 

and the cubic term in the slow roll approximation [27] is 

S3 = J dr/dW[-V:,j^50^(50f50f + 5,J5A,^^50^'9V-2(50f')'9^0f- 
6ijSMN^S<P^S<p{'<Pi' - 6jj6MN^S<l)^d6<l)id6<p(] + perms , (20) 

ori ori 

where the permutations represent swapping the 502 in other possible positions. The 
term containing Vjjk in the action is not neglected by the slow- roll approximation, 
because it may contribute large effect in non-Gaussianity [23j. 

Variation 5S/5{5(j)^2) — gives the field equation. All the surface terms are ne- 
glected, which requires that 5(50() vanishes in the boundary, and the equation of 
motion for simplifies the result further. The final result is 

(50^" + 2m0^' + P50^ 

= {-a'VuK54>i5ct>^) + ^{-25MNdV-^5(t>'^'d6(t>{' + 25MNH'^VHct>{) 

ri 

+V250f V2(50f + (50f' V2(50f ' + 9(50f' 550^' )] • (21) 

On the right hand side of the equation it is the source term. Using the Green's 
function, takes the form 

^</'2(^, x) = I (0e^^- I dr j '^^0^Gr,{r^^ t)5{c^ - k^ - k,)5| , (22) 
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where 

ri rt 

+5M^^i^50f' + 5Miv^2^50f , (23) 
rt ri. 

and {JFi, JF2, ^1, Q2} are some factors in the momentum space. 

T, = -2kl, = Iki . k, - ^ (kfkl + tjki . k,) (24) 

ft = -|k, . k. . = 4 + ^3^^ (fe! + k, . k.) (25) 

Notice that the terms with JF2 and Q2 are symmetric with M, A^, so when we calcu- 
late the three-point function, JF2 and Q2 must be symmetrized over permutations of 
{^1,^2} as in p.8j. On the other hand, J-'i and Qi cannot be symmetrized. 

5 Three-point correlator 

The three-point correlator of a free scalar field vanishes, = 0. The leading 

order of three-point correlator {S(j)S(f)S(l)) is (50i50i502) ~ |(50i<^0i50i * where 
* denotes a convolution. Thus with the value of 6(j){ iHM and 502 ( |22l) . the three- 
point correlator of multiple- field can be calculated. As the argument given in |28j . 
the field equation of multiple- field is in the approximation of slow roll limit, and the 
expansion in powers of slow-roll parameter is not applicable at the end of inflation. 
The reason is that the the subleading term has logarithmic divergences ln|fc?7| = A^, 
and the growth of the e-folding number makes the subleading terms not negligible. 
Here we just calculate the three-point correlator when the modes cross the horizon. 

According to the source term of the field equation, we will show the results of 
three-point correlator from the three parts below. 

5.1 VjjK terms 

In slow roll approximation, we neglect the Vi and Vjj terms, but the Vjjk terms 
could have non-neglect able effect in some situation, and also lead to the logarithmic 
divergence. 
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The three-point correlator from Vjjk terms take the form 

J —CO 

[OkAryekAv) -dkAr)ekAv)\ K{v)OkMdkAr)K{r) + 

K {r)K iv) - K {r)dt, iv)] Ok, {V)0k, {r)ek, (r) | . (26) 

The 9 terms in the brackets come from the Green's function. Outside the brackets, 
the 9 terms which depend on time parameter rj are derived from the free scalar field, 
and the 9 terms depending on r is derived from the source of the field equation. We 
have stated that 9 is independent of the index in multiple-field, so the final result is 
similar to single field for V,i,jk term. With the value of 9 f|T6|) . this part of three-point 
correlator can be written to the leading order of slow roll approximation. 



(50^(ki)50-^(k2)50^(k3)) 3 (27r)35(ki+k2 + k3)-^ 



,IJK 



% Re[-z(l - ihr){l - ik2T){l - ihT)e'^'^] , (27) 

where fcj = fci + ^2 + ^3, * denotes the value at the time rj. Here we take rj to the 
value that the modes cross the horizon. Following the paper [5], we can deform the 
integration variable r to Euclidean time and deal with the divergence of the integral 
properly. Finally, we obtain 



(50^(ki)<50-'(k2)<50^(k3)) D (27r)3<5(ki + k2 + k3; 



UK 
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X 



(-h^ + h]lkk, + ^{^ + ^ + ln\k,v\]j2kf) ' (28) 

where i G {1,2,3}, and 7 ~ 0.577 is the Euler's constant. In multiple-field, there 
also exists the infra-red divergence term form the Vjjk terms. When we take the 
time crossing the horizon, the divergence term is negligible. The classical evolution of 
perturbation afterwards will make the term large, but we can use other formalism to 
deal with the problem, such as 6N formalism, or the separate universe approach 

El. 
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5.2 T terms 

In this section, we discuss the zero-derivative terms in the source, which contain the 
contributions of the T\ and terms. Since there is a delta function in the three- 
point correlator, the sum of momentum ki + k2 + ks is zero, we could write T\ and 
JF2 in other equivalent form, 

Ti{ki,k2;h) = -2kl , (29) 
and JF2 could be symmetrized as 



^2(A:i, k,- k,) = -^ikl + kl) + 1^1^ + ^ . (30) 
The part of three-point correlator come from the J^i and terms is expressed as 
{dcj>\k,)Scj)'{k2)Scj)''{k,)) D t{2n)'5{k, + k2 + ks) r dr airf x ^ x 

[V2e'6'''j^i{k^, k2; k^) + V2e' d'"" Tiik^, k^) + 2V2e^5^-^^2(fci, k2] k^)] 

[V2e'6'''j^i{ki, ks; A^s) + v^e^^d'^J^iik^, k,; k2) + 2v^e'^5^^^2(A;3, ^i; ^^2)] 
[V2e-^5'^J-i(A;2, k^; k^) + V2e^5'-'^i(A:3, ^2; A;i) + 2y2e^5-'^^2(A;2, k^; k^)] 

[OkAr)0kAri) - K{r)ekAv)] OkMOuMKiryeu^ir)^. (31) 

where origins from the commutation relation of creation and annihilation opera- 
tors f|T5l) . and the factor 1/2 ahead of the open brace comes from the definition of 502- 
Different modes of contracting a free scalar field 5(f)i and in the source term make 
JF2 symmetric, and cause the factor 2 in front of J^2- Here the slow-roll parameter 
is defined in (|H]). 

In the three-point correlator, terms are independent of integration variable, so 
the factors of final results /i, /2, and /a are similar to the case of single field, 

(5(/.^(ki)50^(k2)50^(k3)) 3 (27r)35(ki + k2 + ka)^^ x 1 

/i[V2e'(5'^^>i(A;i, k2] k^) + y2e-^5^^^i(A;2, A;i; k^) + 2y2e^5'^J-2(A;i, A;2; ^3)] + 
/2[V2e'5-'^^i(A;i, A^g; ^2) + V2e^5"j-i(fc3, ^2) + 2V2e-'5^^J-2(A;3, fci; A^a)] + 
/3[v^e^5^^^i(A;2, A^s; k^) + y2e^5"^i(A;3, A;2; fci) + 2y2e^5^^^2(A;2, A;3; A;i)]}(32) 
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where 



(kl + k2- 






2kl{kl - Akik^ 


+ kl- 


kl) 


{kl - {k, - 


- k^YY 




2kf{kl + Akiks 


+ kl- 


kl) 


~ {kl - k2- 


ksm 





(33) 



The factor origins from 6 and the scale factor a. 



5.3 Q terms 

Similar to the derivation of T terms, Q terms lead to another part of the three-point 
correlator. We have 

^2 _L ^2 _ 1.2 

If 

and Q2 ~ when symmetrized. The expectation value of Q terms is 



Qi{ki,k2,ks) = ^—^, (34) 



X 



(50^(ki)50^(k2)50^(k3)) 5 z(27r)=^5(ki + k2 + k3) r dTa(T)2 X ^ 
{\^^2e'b''^Ql{ku k2] ks) + V2e''5'''gi{k2, ki; k^)] 

[9k,{r)0kM - 0k,{T)9kM] OkAv)OkM-^o,,{r)^e,,{T) + 

[y/2e'S''''gi{ki, ks; k2) + V2e''S''gi{ks, k^, k2)] 
[OkAr)OM - OkAr)OkM] OkAri)Oks{v)^OkAr)^Ok,{r) + 
[V2e'5'''gi{k2, ks; ki) + V2e'' 5''gi{ks, k2; k^)] 

[OkAr)OkM - K{r)OkM] OkAri)OkMj:^OUr)j:^K{r)}. (35) 

Repeating the progress in the case of T terms, we obtain 

(50^(ki)50^(k2)50^(k3)) D (27r)=^5(ki + k2 + ka)^^ x ^ x 

[gi[V2e'5'''gi{ki, k2; ks) + ^e' 5'''gi{k2, k^; ks)] + 
g2[V2e'5'''gi{ki, ks; k2) + V2e''S''gi{ks, k^; k2)\ + 
9s[V2e'5'''gi{k2, ks; fci) + ^^e"" 5''gi{ks, ^2; fci)]} , (36) 
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where 

9i 

92 



2 



(kj - {k, - ksW 



Then we find the summation of the three-point correlator. For simphcity, the Vjjk 
terms are neglected. 

(50^(kO50''(k2)50^(k3)) 3 (27r)3<5(ki + k2 + ka) X 

^ E ^-/''Aik^k^k,), (38) 

1 i« I perms 

Aik,, k,, k,) = hi - MM±M) _ ^klkl , (39) 

where P,, = H^/Att"^ is the power spectrum of scalar field when modes of scalar field 
perturbation cross the horizon, and A{ki, k2, k^) is the shape factor of the three-point 
correlator. The result is easily reduced to the case of single field, and consistent 

with [SIEB]. 



6 Non-Gaussianity /nl 

The non-Gaussianity of multiple-field inflation is shown in this section. The inflaton 
will decay at the end of inflation, and the flnal observable is the curvature perturbation 
C- To discuss the non-Gaussianity, the power spectrum and bispectrum of curvature 
are given, 

(C(ki)C(k2)) = (27r)35(ki + k2)^Pc(A;i) , (40) 

(C(ki)C(k2)C(ki)) = (27r)35(ki + k2 + k3)Sc(fci,fc2,fc3) . (41) 

The non-Gaussianity is the deviation from the Gaussian statistics in CMB, and use 
the parameter /nl to represent its magnitude, 

C = C. + ^/nl K - (C.')) , (42) 

^ Since we have used the metric of field space J"'^, the resuh is a httle diff'erent from 21] . 
Therefore, /c2 and are symmetric here. If we symmetrize k2 and k^, in eq.(69) of [27j . then we get 
the same resuU. So both results are equivalent. 
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where (g denotes the Gaussian part of (. □ Using this definition, we obtain 

5/NL - ^4;^ • (43) 

Notice that the field equation formahsm which is used to calculate the non- 
Gaussianity is applicable when the modes cross the horizon. In multiple-field infiation, 
there exists entropy perturbation. Thus in order to consider the effects afterward, the 
5N formalism is a good method. On large scale, the value of curvature perturbation 
is the e-folding number from the initial fiat slice at to the final uniform density 
slice at time t, 

C(t, x) ~ (5iV = N{t, t„ x) - N{t, U) , (44) 
where the e-folding number is defined as 

N{t,U) = j^Hdt . (45) 

5N can be expanded by the initial scalar fields, 

5N = Nj5<p^ + ]^Njj5(p^5(l)-^ + ■■■ . (46) 

The power spectrum and bispectrum can be expressed by 5N formalism, 

= 5''NjN^jV, , (47) 

(C(ki)C(k2)C(k3)) = iV,,iV,jiV,^(50^(ki)50^(k2)50^(k3)) + 

]^NjN,jN^kl{H\\^iW{\^2){H'' * 50'')(k3)) + perms , (48) 

where * denotes a convolution and the higher order terms are neglected. The non- 
linear parameter /nl is derived from (H3l) . (H71) and (HHl) . 



The last term on the right hand side is from the curvature evolution on large scale, 
which contributes the local form of non-Gaussianity. The equilateral shape of non- 
Gaussianity from multiple-field is constrained by the tensor-to-scalar ratio r ~ P^/P,^. 



^Here the sign of /nl is consistent with the convention of the CMB experiments, and different 
from the paper written by Maldacena [5j. 
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7 Conclusion 

In this paper, we derive the second-order field equation of multiple-field, and calculate 
the shape of non-Gaussianity for multiple-field infiation with the method of the field 
equation. The shape of Non-Gaussiantiy derives from the three-point correlator, 
which implies the microphysics in the period of infiation. Our result of the three- 
point correlator is consistent with the previous one [27] which uses the method of 
in-in formalism. And it is easy to extend the Bunch-Davies vacuum to the a vacuum 
in the field equation formalism which shows the trans-Plankian physics ^3] from 
non-Gaussianity. 

The field equation formalism is applicable when we know the equation of motion, 
even if the action is not given. Meanwhile, we should notice that the formalism is 
used when the modes of scalar field perturbation crossing the horizon. After crossing 
the horizon, the quantum fiuctuations become classical due to decoherence [36]. Then 
the classical evolution of curvature perturbation on large scale could be solved by the 
5iV formalism. Finally, we could get the non-Gaussianity observed in the CMB. 
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